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Abstract 

We propose a scheme to engineer an effective spin Hamiltonian starting from a system of elec- 
trons confined in micro- Penning traps. By means of appropriate sequences of electromagnetic 
pulses, alternated to periods of free evolution, we control the shape and strength of the spin-spin 
interaction. Moreover, we can modify the effective magnetic field experienced by the particle spin. 
This procedure enables us to reproduce notable quantum spin systems, such as Ising and XY 
models. Thanks to its scalability, our scheme can be applied to a fairly large number of trapped 
particles within the reach of near future technology. 

PACS numbers: 03.65.-w, 03.67.Ac, 03.67.Lx, 75.10.Jm 



I. INTRODUCTION 



Single electrons confined in Penning traps may represent a valid, experimentally viable 
system for the implementation of a quantum processor . Our proposals have been 

encouraged by the astonishing results obtained in high precision experiments with a single 
electron and by the advances in trapping technology, from micro-traps [9| to 



scalable open planar Penning traps 



10 



11] . In this spirit, it has also been put forward how 



to realize a quantum information channel, based on interacting spin chains, by means of 
trapped electrons [ijl, Q . 

In this paper we focus on a linear array of electrons, each one confined in a micro- Penning 
trap. Our aim is to prove that, from the same physical system, we can derive a variety of 
interacting spin models. In particular, we show how to design and control the relevant 
terms in the effective spin Hamiltonian. As a result a system of trapped electrons can be 
exploited to study the dynamics of a wide range of quantum spin models. We recall that 
these models are very important for the understanding of the rich phenomenology observed 
in several quantum many-body systems, such as quantum magnets and high temperature 
superconductors. Moreover quantum spin systems are able to exhibit quantum phase tran- 
sitions lij]. To this end, it is critical to control and vary system parameters like the applied 
magnetic field and the spin-spin coupling strength. Our method to shape the effective spin- 
spin interaction employs sequences of electromagnetic pulses alternated to periods of free 
evolution. This technique is similar to the refocusing schemes used in nuclear magnetic 
resonance (NMR) experiments and relies, from the theoretical point of view, on the average 



Hamiltonian theory 



1 51 ] . We point out that a system of trapped electrons presents several 



advantages over NMR implementations. The most important ones are scalability and the 
possibility to independently adjust the values of relevant quantities, like the spin precession 
frequencies and the spin-spin coupling. Indeed their values depend on external parameters 



such as the magnetic fie 



the trap electrodes 



12j. 



d gradient, the inter-particle distance and the voltage applied to 



Also other systems, like linear or planar arrays of trapped ions, enjoy some of these 
properties and, therefore, it has been proposed to use them as a quantum simulator for 
interacting spin chains \m, \v\, 18]. However working with trapped electrons we naturally 
have a system of spin one-half particles, without the need for artificially creating an effec- 
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tive two-level system. Another difference between trapped ions and electrons relies in the 
typical resonance frequencies. Ions are controlled by means of a sophisticated laser setup, 
while trapped electrons are manipulated by microwave or radio-frequency fields. In this 
respect trapped electrons can benefit from the same technology already developed for NMR 
spectroscopy. 

In th is p aper we consider both the case of electrons with the same spin precession fre- 
quency [l2 1 and the case of electrons with different spin precession frequencies jsj. The 
frequency addressability, which is necessary to manipulate specific particles in the array, is 
obtained with the insertion of a magnetic field gradient. However this condition is required 
only to modify the interaction range and topology. In the case of electrons with the same 
spin precession frequency, we prove that, by flipping the spin state twice, we can effectively 
reduce or even cancel the spin dynamics due to the external uniform magnetic field. This 
way the effective spin system is subjected to a weaker magnetic field, whose intensity can go 
down to zero, without affecting the overall trap stability. The same resonant electromagnetic 
field, used to flip the spin, is able to produce coherent superpositions of the two spin states 
by adjusting its phase and duration. These operations are the building blocks of specific 
pulse sequences that allow to engineer the effective spin Hamiltonian. By iterating such 
pulse sequences, we obtain various interesting spin Hamiltonians such as the Ising model or 
the XY model. In addition, if we want to customize the interaction range and coordination 
number, we should apply similar sequences of pulses to selected subsets of spins in the ar- 
ray. The resulting spin system can exhibit a nearest neighbor (NN) as well as a long range 
interaction. The number of pulses in each sequence is relatively small and, most notably, 
does not depend on the number of spins in the array, thus making our procedure scalable. 

The paper is organized as follows. In Sec. [Til we briefly present the system of trapped 
electrons and review the derivation of the effective spin-spin interaction. In Sec. II III we 
describe how to prepare and manipulate, with an additional oscillating magnetic field, the 
spin state of each electron in the array. In Sec. IIVI we show how to engineer the spin 
Hamiltonian by applying appropriate electromagnetic pulse sequences, that allow to control 
the strength and the range of the interaction. The capability of our technique to reproduce 
a given Hamiltonian is analyzed in Sec. |V] Finally in Sec. I VII we summarize our results and 
discuss future perspectives. The more technical details, concerning the design of the pulse 
sequences and the estimate of fidelity, are reported, respectively, in Appendices |A] and [B] 
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II. ARRAY OF TRAPPED ELECTRONS 



Let us consider a system of N electrons confined in an array of micro-Penning traps in 
the presence of linear magnetic gradients. The Hamiltonian of the system can be written as 

N N 

H = T,Hr + T,H^ (1) 

i=l i>j 



where 

jjNC _ (Pi 

2m P 1 4m 



e 



represents the single electron dynamics inside a trap and 



H u = 7 — r~ r (3) 

,J 47re |ri-rj| 

describes the Coulomb interaction between electrons i and j. In Eqs. (J2]) and m e , e, 
g, and er^ are, respectively, the electron mass, charge, gyromagnetic factor, and Pauli spin 
operators. We assume that the micro-traps are aligned along the x axis and that x^q is the 
position of the center of the i-th trap. The electrostatic potential 

fe,*.)^.^^* (4) 

is the usual quadrupole potential of a Penning trap, where Vq is the applied potential dif- 
erence between the trap electrodes and i is a characteristic trap length. The magnetic field 

Bi = --[{xi - ar i>0 )i + yd] + (B 0:i + bz^k (5) 

is the sum of the trapping magnetic field B 0i h, providing the radial confinement, with a 
local linear magnetic gradient b around the i-th trap. The associated vector potential 

A; = ^(B 0>i + bzi)[-yii + (a* - x i)0 )i] (6) 

preserves the cylindrical symmetry of the unperturbed trapping field. 
Following the approach described in 

ay 

. Il2l | the Hamiltonian, Eq. (j2J), of a single electron 

can be written as 



+ -e%u z (a Z)i + 4 :i ) a- - ^ehuj z J^- (a\ +) a c ,i + a\ (7) 



where the annihilation operators a mj i, a c> i, a z ^ refer, respectively, to the magnetron, cyclotron 
and axial oscillators of the i-th electron and = (erf ±zcrf)/2. The frequencies of the dif- 



ferent electron motions are ~ uj 2 z /{2uj c ^) 1 uj Cji ~ (\e\B 0ti /m e ) —<jj m ^ uo z = J2eV^/ (m e £ 2 ) 
and u S) i = g\e\Bo ) i/(2m e ). The Hamiltonian ([7]) has been obtained under the assumptions 
u m ,i Cw z < uj C) i and 6| ^ | / -Bo,?. *C 1- We also assume that the cyclotron motion is in the 
ground state and the amplitude of the magnetron motion is sufficiently small (axialization) 



19(| . The dimensionless parameter 



efr / ft ,„ s 
m e o; 2 V 2m e oj z 



represents the coupling, due to the magnetic gradient, between internal and external degrees 
of freedom of the particle. 

Similarly, if the oscillation amplitude of the electrons is much smaller than the inter-trap 
distance, the part of the Hamiltonian describing the Coulomb interaction can be written as 

GT 



H ?j - h kj( a *,i + a U)( a zJ + a lj) - ( a ^ a h + a U a ^j) > ( 9 ) 



where = e 2 / (87ieom e u z d^j) with di t j being the distance between the i-th and j-th particle. 
Now we apply to the system Hamiltonian the unitary transformation [20] 

JV 
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8=1 * . 



of (4,i - + —J— (af } ali - crJ +) a C;i ) , 



(10) 



with u at i = u s> i — o; C) j. This transformation formally removes, to the first order in e, the 
interaction between the internal and the external degrees of freedom in Hamiltonian (J7J) 
and, at the same time, introduces a coupling between the spin motions of different electrons. 
Consequently the spin part of the system Hamiltonian can be recast as [12] 

N h fi N 

# S -E • ttE (2^<o| - j;y;-aj - , (11) 



1=1 Z Z i>j 



where 



f g\ 2 f 9\ 2 ftc 4 ^ 2 

J ^ = I2J Ci>j£ = {2) 167ie mi Utdf/ (12) 
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The effective spin Hamiltonian (fTTT) exhibits a long range interaction between all the particles 
in the chain. The coupling strength decreases with the third power of the distance between 
particles, i.e. with a dipole-like behavior. Moreover, Jfj and depend, respectively, on the 
axial frequency and the cyclotron and anomaly frequencies. Since the trapping frequencies 
form a well defined hierarchy, the coupling in the longitudinal and transverse direction can 
be utterly different. For example, for typical experimental values of the cyclotron and axial 
frequencies, such as u c /2n ~ 100 GHz and u z /2tt ~ 100 MHz, the ratio J^J/Jf^ is less 
than 10~ 6 . Therefore, for practical purposes J^J is often negligible with respect to J?-. In 
particular, this is true when the difference between the spin frequencies of different particles 
is much larger than their xy spin-spin coupling strength. In this case the spin Hamiltonian 
reduces to 

N n N 
ff.^E^+AEW ( 14 ) 

i=l Z i>j 

In Hamiltonian Eq. fTH|) we have used the rotating wave approximation (RWA) to neglect 
the interactions between spins along the x and y directions, since they give rapidly rotating 
terms. The Hamiltonian (HM is, therefore, simi 
molecules used to perform NMR experiments 
frequency differentiation and the spin-spin couplings are determined by the chemical nature 
of the molecules, whereas in our system they depend on the value of the applied fields, that 
are under control of the experimenter. 



ar to the nuclear spin Hamiltonian of the 

n 

15| . However in NMR systems the spin 



III. SPIN STATE MANIPULATION 



In this section we describe how to prepare and manipulate the spin state with an external 
oscillatory field. Let us consider a magnetic field h p (t) oscillating in the xy plane with 
frequency u and phase 6 such that 

b p (t) =b p [icos(ujt + 6) +j sm(ut + 0)]. (15) 

If we add this field to the system, the spin Hamiltonian, Eq. (iTTj) . becomes (here and in the 
rest of the paper we set ft — 1) 

1 n N 
H ^ 5 E "s,i°j + | EK (+) e" lM+9) + c7^e^ t+ % (16) 

1 3=1 Z 3=1 
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with x = g\ e \b P /{2m e ). In deriving the Hamiltonian (flEI) . we assumed that the interaction 
between the electrons and the oscillating magnetic field is much stronger than the spin-spin 
coupling. Hence, the terms in Eq. (TTTT) proportional to Jfj and J^J can be neglected. In the 
case of a system with spin frequency differentiation the field ( I15p . applied for an appropriate 
time t with frequency to = uj s j, affects only the spin states of the resonant j-th electron 

I T>i - e^-/^ cos (^) I !>i - z e ^-'/ 2 )*+*> sin (^) | |),, (18) 

Without spin frequency differentiation, the single qubit addressing with microwave radiation 
is, of course, no longer possible. Therefore when all the spins have the same precession 
frequency u s , a single resonant pulse suffices to produce the evolution of Eqs. (fTTj) and ffTB"]) 
for each particle in the array. From Eqs. ( fT7l) and ( fl8l) we see that by changing duration 
and phase of the applied pulse we can prepare and manipulate at will the spin states of the 
trapped electrons. In particular, if we apply a pulse for a time i — tt/x with 9 = 0, we can 
flip the spin state of each particle 

I i)j -> -te-^^l T)„ (19) 
| t)j -> -ie*^/ 2 ^! I),, (20) 

We define this transformation as 

TV 

F = (g){-i[(jj +) e- i( ^^ /2) * + a^e^l^ 1 }}. (21) 
i=i 

It is not difficult to verify that the inverse transformation F~ l is obtained with a pulse of 
the same duration i but with phase it. 

If we move to the interaction picture (IP) with respect to the Hamiltonian Y,iLi( UJ Si i/2)<7i , 
the system evolution is given by Eqs. f|T7|) and ffTSl with u s j = 0. Consequently, the spin 
flip operation, Eqs. (TT9l and ( 1201) . turns into 

I 1), - -<| t>i, (22) 
I T)i ^ -<l Di- (23) 

The above transformations correspond to the application of the operator —icTj- In a similar 
way a pulse applied for the time t with 6 = n/2 produces a transformation corresponding 
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to the application of Furthermore, always working in IP, if the pulse is applied for a 

time i/2, we can obtain the pseudo-Hadamard operations 

N (1 - ia x ) 

G x = (g) ^ for 9 = 0, (24) 



V2 



y (1 + ia x ) 



V2 



G y = ®^^ for0 = -, (26) 



>/2 2' 

for ^-f. (27) 

The coherent superposition of the spin states | |), | |) for each particle can be achieved with 
a single multi-frequency pulse. Hence, an appropriate choice of the frequency, duration and 
phase of the pulses allows for performing, apart from irrelevant phase factors, single qubit 
operations on each spin of the array. 



IV. ENGINEERING THE SPIN HAMILTONIAN 

In this section we show that, by using the additional magnetic field ( jl~5{) . we can also 
adjust and control the form of the effective spin Hamiltonian, starting from the models 
given by Eqs. (TTTT) and (fl4l) . This is achieved by applying to the system specific sequences 
of pulses alternated to periods of free evolution. Our approach is inspired to the refocusing 
schemes used in NMR experiments [l5| . Similarly to this technique a key point is the choice 
of the different time scales. Spin operations, operated by means of pulses, should be virtually 
instantaneous with respect to the free evolution of the system. Therefore, the pulse duration 
should be much shorter than the free evolution time. 



A. Tuning of the effective magnetic field 



The spin Hamiltonian, Eq. ( JTlj) . in the case of spins with the same precession frequency 
can be recast as 

H s ~ H + H c , (28) 
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where 

N 



Ho - E (29) 
i=i z 



1 

tf c = 2 E - J >J a ^ - J » v d ■ ( 3 °) 

In the following we shall prove that, by sending resonant pulses of the kind of Eq. (|T5l) . 
it is possible to reduce or even cancel the effects on the spin dynamics of the Hamiltonian 
term Hq. This result corresponds to an effective modulation of the external magnetic field, 
without affecting the trapping stability of the whole set up. 

In particular, by applying a sequence consisting of a pulse producing the spin flip trans- 
formation F, Eq. (T2T]), followed by a period of free evolution t and by a pulse producing the 
inverse transformation -F -1 , we can change the sign of the Hamiltonian term H 

F -l e -iH s t F = exp [_^_ Fo + £Q t ]_ ( 31 ) 

To prove Eq. fl31|) we use the identity 

F -i e - iHa t F = expf-^F-^o-F + F~ x H c F)t\. (32) 

Now we have 

N /, , \ 

F- l H F = f^a^e-^/^+a^e^/^] [^e^'^ + a^e^'^] 

j= i V 2 / 

= -Ey°? = -ffo. (33) 
i=i 

The identity F~ l H c F = H c follows from the commutation relation [H c , F] = 0, which can be 
verified with some algebra. Moreover we observe that [Hq,H c ] = 0, because the interaction 
Hamiltonian preserves the total magnetization Y^=i^t- From this last consideration and 
from Eq. fl3Tl) we find 

F -x e -M.t* Fer iM.ti = exp [_^ eff (t 1 + f2 )] > (34) 

with 

H eS = t --^H + H c . (35) 

The left hand side of relation (13^1) represents a sequence consisting of a period t\ of free 
evolution, a pulse producing the transformation F, a period t 2 of free evolution and a pulse 



producing From the right hand side of Eq. fl34l) . we see that this sequence is equivalent 
to the system evolution for the total time t\ + t 2 according to the Hamiltonian if e g. Hence, 
we can obtain an effective reduction, by a factor (ti — t 2 ) / {ti + £2), of the Hamiltonian term 
H . This result can be viewed as a decrease of the magnitude of the uniform magnetic field 
as far as the electron spin dynamics is concerned. Notice that for t\ = t 2 we can completely 
suppress the dynamical effects due to the term H . 



B. Design and control of the spin-spin coupling 

Let us now consider a system with spin frequency differentiation. If we add another field 
consisting of a superposition of terms resonant with the spin frequencies 

N 

b s (t) = b s [i cos(u s , k t) + j sin(u Syk t)}, (36) 
fc=i 

the spin Hamiltonian of Eq. (Tl4"l) becomes in IP with respect to X^i(^s,i/2)crf 

jjiP ~H Z + H bs , (37) 

with 

JV 

#* = £^M> (38) 

i>j 

N N 

H bs = vT,{^ ) +^ ) )=^h (39) 

i=l i=l 

where rj = g\e\b s /(4m e ). Hence, the application of the oscillating field (I3"B"|) gives rise to an 
effective static transverse magnetic field, whose strength can be controlled and modified, 
since it depends on the field amplitude b s . This tool may turn out useful in reproducing 
quantum models like Ising system of spins. In this case the parameter rj should be comparable 
to the coupling strength J[j between the spins. 

Moreover, we can engineer the spin-spin coupling, that is introduce an effective spin-spin 
interaction along the x and y axes. This is achieved by means of sequences of pulses, of the 
kind given in Eqs. (1251) . fl2B"j) . and (j27|) . affecting all the spins in the array. Indeed, it 

can be easily proved that 

G x e~ iHH Gl = e - lG * HZGU = exp l-i J^ajtxft ] , (40) 
G y e- iHH Gl = = exp ( -i £ J* k (T*a%t ] . (41) 
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Hence a sequence of two specific pulses, alternated to a period of free evolution under the 
Hamiltonian H z , effectively modifies the direction of the spin-spin coupling. Now, if we 
combine the three operations (1551) . (jlQ|) . and (T4"Tj) we have that for ti,t 2 ,t 3 <C n / J[ i+1 , n / 1] 

e-^ HZ+Hbs ^(G y e~ tH ^Gl)(G x e- iHHs Gl) ~ e -^(*i+*»+*s), (42) 



with 



V V AT V 

^eff = ri77 J2 °i + n E J 4>M + r 3 E J i.j CT : (T j- ( 43 ) 

i=l 4>j i>jr 4>j 



where r, = t»/(ti + £2 + £3)- In deriving relation f|42l) we used the approximate identity 21| 

e -jAltl e -iA2<2 g-i^nin r^, g — i (Tl Al +T2 A2 + . . .+r n A n )t (44) 

with £ = Y^i=\ti an d Tj = tj/t, which is valid, to first order in t, for ti much shorter 
than the typical time scale of the dynamics due to the Hamiltonian Aj. However, more 
elaborate sequences of pulses (see Appendix [B]) give approximations to higher orders in t 
22I . A recursive application of the sequence (142]) determines an effective evolution under the 
Hamiltonian H e R. We point out that, in this Hamiltonian, Eq. (j43l) . we can independently 
control and change the values of the parameters Tj's and rj, since they depend, respectively, 
on the free evolution times tj's and on the pulse amplitude b s . Consequently we obtain an 
Hamiltonian H e ^ with a variable relative strength of the spin-spin coupling in the x, y, and 
z directions and a tunable transverse magnetic field. Notice that we can also set 7* = 
for any desired % or rj = 0. This is achieved by simply choosing ti = or switching off the 
external field h s (t). In this way various interesting quantum spin models can be derived 
from Hamiltonian (j42p . For example for r 2 = r 3 = we obtain the Ising model, whereas for 
Ti = (t 2 = 0) we obtain the XY model in its usual (rotated) basis. 

C. How to modify the interaction range and topology 

The spin Hamiltonian (1381) can be written as 

N 

E =H Z 1+ HI + ..., (45) 

i>j 

where 

N-n 

K=H Jti+n^hn (46) 
4=1 



ff 



represents the coupling between the n-th nearest neighbor spins. In our system the inter- 
action between spins has a dipole-like nature, i.e. it decreases with the third power of the 
inter-particle distance. Consequently only the first few terms H* at the right hand side of 
Eq. (jl"51) play a significant role. In the following we are going to outline a procedure to in- 
dependently control and modify, in a relatively simple way, the strength and the sign of the 
relevant terms. In other words we can design the interaction topology by enhancing or 
suppressing the coupling between the n-th nearest neighbors. This is achieved by iteratively 
applying to the system appropriate sequences consisting of pulses alternated to periods of 
free evolution. In our scheme each pulse affects simultaneously a specific subset of spins in 
the array. Notice that if a particular pulse sequence S modifies the spin Hamiltonian H z , we 
can extend the same kind of coupling to the other directions by simply performing, according 
to Eqs. (j4"U|) and (jHJ), the sequences G y SG^ y and G X SG\.. Therefore, we restrict ourselves 
to the transformations affecting the spin-spin coupling along the z direction. Furthermore, 
we are going to prove that the number of pulses in each sequence does not depend on the 
number of spins in the array, thus making our technique scalable with the system size. 

As an example, we describe how to suppress the second nearest neighbor interaction by 
making use of three different transformations, defined as a x , a® and <fL. Each of them 
can be performed by a single multi-frequency pulse. The transformation a x consists in the 
simultaneous application of o x to all the spins in the odd sites of the array. The transfor- 
mations o x Cl and a x 2 flip, instead, alternated couples of neighboring spins. In particular, a x ± 
affects the spin couples {1, 2}, {5, 6}, . . ., whereas a x 2 affects the couples {2, 3}, {6, 7}, . . .. 
We prove in Appendix |A] that the sequence 

a x 2 e~ iHZ t a x e~ iHZ I < e~ iHH (47) 

corresponds to the system evolution for a time t under the effective Hamiltonian H c g ~ 
HI+H3, where the coupling of each spin with its second nearest neighbors has been removed. 
Indeed, since the term if| is small, the above sequence well approximates an effective NN 
Hamiltonian. The transformations a x , a x Cx and a x 2 are the building blocks to construct other 
sequences which realize different kinds of spin Hamiltonians. For example, as described in 
Appendix|Al we can easily invert the sign of H z , thus switching from an anti- ferromagnetic to 
a ferromagnetic interaction, or make the coupling strength between first and second nearest 
neighbors equal. This last case corresponds to an effective change in the array topology, 
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FIG. 1: Schematic drawing of a linear chain of N spins (upper part). When the first and second 
nearest neighbor coupling strengths are made equal, the linear chain becomes equivalent to a planar 
array (lower part). 

since the number of nearest neighbors passes from two (linear chain) to four (see Fig. [1]). 

In order to affect the coupling between neighboring spins of order higher than three, we 
should apply simultaneously the transformation a x to selected subsets of three or more spins 
along the chain. For example, as we show in Appendix [A], with a seven pulse sequence we 
can suppress the interaction between a spin and all its neighbors from the second up to the 
sixth nearest neighbors. In such a way we improve our approximation of a NN interacting 
spin chain. 



V. FIDELITY 



In this section we discuss the performances of the scheme, based on the detailed analysis 
reported in Appendix [Bl We emphasize that our treatment focuses on the limitations due to 
the mapping of the system of trapped electrons into the desired target system of interacting 
spins. Hence, most experimental imperfections are not considered here. 

As described in the previous section, to derive the effective spin Hamiltonian we make 



use of the approximate identity (Hl|) or of more sophisticated approximations [22]. As a 
consequence we introduce an error [21] 



E=\\U - U'\\ = max m . m[=1 \(U - U')\if))\, (48) 

which measures the distance between the desired evolution U and the approximated evolu- 
tion U'. For instance, in our case the target unitary operator U = exp(—iH e gt), with the 
effective spin Hamiltonian H c r of the kind of Eq. (|43|) . is approximated to the fourth order 
in t by the sequence S, Eq. (1B2|) [22j. By using some algebra (see Appendix [B]) , we can 

13 



bound the error associated to the application of a single sequence S with 

S s < (J z t) 5 f(N), (49) 

where J z = Jf i+1 is the NN coupling strength and f(N) is, in good approximation, an 

increasing linear function of the number N of electrons in the array. The exact form of f(N) 

depends on the specific target spin Hamiltonian. From this result, we see that the error is 

small whenever the time evolution is much shorter than the flipping time, i.e. t <C n/J z . We 

prove in Appendix [B] that, if we iterate m times the sequence S the total error is £ < mSs- 

Therefore, if we want to simulate the system evolution for a given time T = mt, to keep the 

accuracy high we should apply the same sequence S m times 

„ (J z Tff(N) , . 

6 < - — : ■ 50 

m 4 

For a given simulation time T and coupling strength J z , the error £ decreases with the 
number of iterations m and, therefore, with the total number of pulses. 

In appendix [B] we provide the explicit expression of f(N) for the XY and NN Ising 
models. Consequently we are able to estimate the upper bound of £ in both cases. In our 
analysis, we also take into account the error introduced by the derivation of the effective 
spin-spin coupling 12) . In particular, the error £ c due to the canonical transformation (flQj) 
satisfies the relation 



£ c <N(k + -je 2 , (51) 

where k is the mean axial oscillator excitation number. We consider an array of 50 electrons 
with inter-particle distance of 100 fim, u s /2tt = 100 GHz, oj z /2h = 160 MHz, and a magnetic 
gradient b ~ 200 T/m. With these parameters we obtain, according to Eq. ffl2l . a NN 
coupling constant J z = 10 Hz. We also assume that the spin frequencies of neighboring 
electrons differ of about 2 MHz, each pulse has a duration of the order of /xs and the axial 
motion is cooled to the ground state. In order to simulate the XY model (NN Ising model) 
for a time T = 1 s with fidelity of 99%, we need to iterate the specific sequence S about 
100 (50) times. In particular the simulation of the XY model requires about 3000 pulses, 
whereas the Ising model with NN coupling requires about 2000 pulses. Notice that the Ising 
model with dipole-like coupling requires no pulse sequence, since it is obtained, directly, by 
applying the field ( |36l) . Therefore, in this case we only take into account, as a source of 
error, the thermal excitation of the axial oscillator which, according to Eq. (15T1) . is of the 
order of 10~ 3 . 
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VI. CONCLUSIONS 



In this paper we have proposed a scalable technique for easily controlling and adjusting 
the effective Hamiltonian of a system of interacting spins. The underlying physical system 
consists of an array of trapped electrons in micro-Penning traps. The electron spin is pre- 
pared and manipulated with an external resonant magnetic field. These spin operations, 
applied to all the particles or subsets of them, are alternated to periods of free evolution in a 
fashion similar to NMR refocusing schemes. To selectively address the electrons in the array, 
it is necessary to introduce a detuning between the characteristic spin frequencies by means 
of a magnetic gradient. In particular, we have shown that, in the case of a system without 
spin frequency differentiation, a two pulse sequence permits to reduce or even cancel the 
effect on the spin dynamics of the uniform magnetic field, without affecting the overall trap 
stability. This is potentially useful for the observation of quantum phase transitions 14|, 
where it is important to modulate the ratio between the external magnetic field and the spin- 
spin coupling strength. In the case of a system with spin frequency differentiation, we have 
proved that with a repeated application of appropriate pulse sequences we can modify and 
control the interaction terms in the effective spin Hamiltonian. As a result a wide range of 
spin Hamiltonians can be obtained, such as the Ising model and the XY models. Moreover, 
specific pulse sequences allow to control the sign and strength of the coupling between the 
k-th nearest neighbors for any significant value of k (first, second, . . . , nearest neighbors). 
As an example, we provide a prescription to obtain an Hamiltonian with substantially only 
NN coupling starting from a dipole-like interaction. In our scheme the number of pulses in 
each sequence is relatively small and does not depend on the number of spins in the array. 
We derive an analytical formula to estimate the fidelity of our method for simulating the 
effective spin Hamiltonian, as a function of the coupling strength, the simulation time and 
the number of particles. Our estimates show that it is feasible to simulate the Ising, with NN 
coupling, and the XY model with fidelity of 99% for a system of 50 electrons with a coupling 
strength J z = 10 Hz. Of course, the evaluation of the performances of a real experiment 
would require a closer analysis of all the possible sources of errors and decoherence. This is, 
however, beyond the scope of the present work. 
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APPENDIX A 



In this appendix, we are going to prove that spin flip operations, applied to subsets of 
particles in the array, result in an effective sign change in the interaction between neighbors 
of arbitrary order. The starting point is represented by the relation 



which reverts the sign of the z-th spin operator. We define the following operators 

N/2 

°2i-l> 



(To 



i=l 



a 



Vitf+i with c fc = {k, k + 4, k + 8, . . .} for k = 1, 2, 



iec fe 



a. 



a 



•'< +1 af +2 with % = {k, k + 6, k + 12, . . .} for k = 1, 2, 3, 



(Al) 

(A2) 
(A3) 
(A4) 



^X +1 af +2 af +3 with = {k, k + 8, k + 16,...} for k = 1,2, 3, 4, (A5) 



i&Qk 



that affect simultaneously different subsets of spins. Moreover, we observe that 



K) 2 = «) 2 = K) 2 = (^J 2 



for any possible value of k, so that we can use the identity 

Ae B A = exp(ABA), 

which holds true for any pair of operators A and B, whenever A 2 = 1. 
From the relation flAll) and the definition of a x Q it follows that 



(A6) 



(A7) 



-l)*+w 



2 

i' 



(A8) 



which amounts to a sign change in the interaction between spins with different parity. Con- 
sequently, given the Hamiltonian H z of Eq. (138]) . the transformation a x H z a x inverts the 
coupling between neighbors of odd orders 

N 

(A9) 



a x e -iHH a x = e -iv* H*v* t = exp 



-*E(-i) fc ^ 

k=l 



This property allows us to make equal in strength the coupling between first and second 
nearest neighbors 

N-2 



r —iH z -t r —iH z t 

a n e s ale ~ 



exp 



-if °j( a hi + °"J+2)* 



(A10) 
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with J' = (2/9)J& +1 . 

With a three-pulse sequence 



T — ?' H Z — T —1 H Z — T 

a c 2 e 2<J o e 2<J c 1 = exp 



N/2 

fc=l 



(All) 



we remove the coupling between odd order neighbors and alternatively change the sign of 
the coupling in H z between even order neighbors. To prove Eq. (lAllj) we use the identity 
of = af,£7? and the commutation relation \a s LH z (fL. a x H z a®] = in order to obtain 

\j C2 ^-1 L ^-2 ^2 ^1 - 

a" e - iH **<7*e- iH **<7? = <f* e ~ iH ^a x a x e~ iH ^a x = e -*(^^+«^ff"«^)|. (A12) 

C2 (J t-2 C2 t-1 t-1 * * 

The transformation a x k H z a x k selectively changes the sign in H z to the operators a z and <rj +1 
with j e Cfe, according to Eq. flA3|) . Consequently we have 



N-l N-2 N-3 

<# 2 < = E - E ^+2« + 2 + E (-iyj?,i + 3°M + 3 + ■ • (Ais) 

j=l j=l i=l 

Af-1 Af-2 AT-3 

= E (-ir-WK+i - E + E (-i) i+1 ^K +3 + • • ( A14 ) 

Hence, to demonstrate Eq. (jAUO . we sum Eq. f1A13[) and Eq. (IA14I) obtaining 

N/2 

a x C2 H z a x C2 + a x c H z a x Ci = -2H Z + 2H Z + . . . = 2 J2(-l) k H* k (A15) 

fe=i 

Notice that in the sum the coupling between nearest neighbors of odd orders cancels out. 

By combining the sequences (1A9j) and OAlip . we can invert the sign of the coupling up 
to third nearest neighbors 

(^e^VA'J^Al - exp [-i{-Hl - H* - H z )t}, (A16) 

thus turning a ferromagnetic interaction into an anti-ferromagnetic one and viceversa. An- 
other consequence of Eq. (lAllj) is 



a x C2 e- iHZt ia x e- lHzt ia x Ci e- iHH ~ exp [-«(#* + Hi)t], (A17) 

where the coupling between second nearest neighbors has been removed. 

The approach described so far can be extended in order to cancel coupling terms of higher 
order. For example, to remove both the second and third nearest neighbor couplings in H z 
we make use of the transformations defined in Eq. (IA4I) . They simultaneously affect sets of 
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three nearest neighbors in alternate succession. With arguments similar to those used for 
verifying Eq. flAllj) . we can demonstrate the following identity 



Ke-^'iofJ^e^iff^Ke^'sofJ = e -i(H;-H S -3H§-HS+...yt/\ ( A 18) 
By combining Eq. flAllj) and Eq. (IA18I) we prove that the sequence 



(■■> 



e i(J o e 3(7 ciH°T 3 e 30 T 3 H°T 2 e 3ff r 2 )i <7 r 1 e 3(T T 1 > e ( A19 ) 



corresponds to the evolution for a time (4/3)t under the Hamiltonian (H* + HI + ...), where 
the coupling between second and third nearest neighbors has been removed. The implemen- 
tation of this sequence requires six pulses, since each couple of consecutive transformations 
in Eq. (1A19j) is equivalent to a single transformation affecting simultaneously a specific 
subset of spins in the array. 

It is worth to point out that with a seven-pulse sequence we can approximate the NN 
model in a very accurate way, i.e. we can suppress the interaction between a spin and all 
its neighbors from the second up to the sixth nearest neighbors. This is achieved by using 
the four transformations defined in Eq. (1A5j) . that simultaneously affect alternated sets of 
four nearest neighbors. It can be proved that 

(A20) 

From Eq. flAllj) and Eq. (IA20I) we have that the sequence 



i^e— *<)K Q4 e- 1 " ^J^e^ ^){p^e^ lo&^e ~ >a" Ql )e 

(A21) 

corresponds to the evolution for a time 2t under the Hamiltonian (Hf + if| + •••)> where the 
coupling between nearest neighbors from the second up to the sixth order has been removed. 
The implementation of the above sequence requires just seven pulses, since each couple of 
consecutive transformations in Eq. (1A21I) is equivalent to a single transformation affecting 
simultaneously a specific subset of spins in the array. 
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APPENDIX B 



The sequence of unitary operators to the left hand side of relation (|44|) approximates the 
evolution under the target Hamiltonian 

n 

E = Y j r i A l (Bl) 
to first order in t. However, more elaborate combinations of unitary operators provide better 



approximations. For example the sequence 22] 



C — CCCC QQQriririririririri CCC (~RO\ 

D = OiDiDiD_2*-'1*-'1*J1*-'1*J1*-'1*J1*-'1*-'1*J1*-'— 2*J1*-'1>-'1 \ ) 

with 

S k = e-'T^e-^* 2 ^ 2 . . . e-^ tnAn (B3) 

and 

S k = e ~ l & nAn . . . e -^ A * e -^ A \ (B4) 

approximates the unitary operator e~ lHt , with t = Y0£=i U and = ti/t, to the fourth order 
in t. 

The error introduced by approximating the unitary operator U with the unitary operator 



U' can be measured by the quantity [21] 



S=\\U- U'\\ = max m . m[=1 \(U - U')\1>)\. (B5) 

Now if we want to approximate the evolution under the Hamiltonian Eq. (IBlj) for a time 
T = mt we can apply m times the sequence (IB2j) . From the generalization of the inequality 



\\AB - CD\\ < \\A — C\\ + \\B - D\\, (B6) 

verified for any unitary operator A, B, C, D, we have that the error £ of our approximation 
satisfies the inequality £ < Yh=i where £i is the error introduced by the z-th application 
of the sequence flB2p . More specifically we evaluate the error 

£ x = We-iT.t^y - S\\, (B7) 

when the operators Aj's are typically of the kind Y,i>j Ji,j a i a> ji with k = x, y, z. We explicitly 
expand each operator to the right hand side and find that their difference is proportional to 
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t 5 , because the sequence S approximates the desired unitary evolution exp(— i J27=i T iAit) to 
the fourth order in t. Moreover, we make use of the inequality 

\\aA + (3B\\ < \a\\\A\\ + \P\\\B\\, (B8) 

which holds true for any pair of operators A and B and complex numbers a, (3. Finally, we 
observe that ||C|| = 1 if C is any product of Pauli operators of. This approach lead us to 
the following estimate for the error, defined in Eq. (1B7j) . 

Si < (J z t) 5 f(N), (B9) 

where J z = J z i+1 is the nearest neighbor coupling strength and f(N) is, in good approxi- 
mation, an increasing linear function of the number iV of electrons in the array, depending 
on the specific form of the spin Hamiltonian. 

When we apply m times the sequence S, from the previous discussion it follows that 
£ < m(J z t) 5 f(N). Now if we indicate with T = mt the total simulation time, we obtain 

(J'TfHN) 

~ m A y ' 

or, equivalently, 



m < fJ^m. (BU) 

Equation ( 1B11I) gives an upper bound to the number of iterations required to mimic the 
desired evolution with an error £ . For example, to approximate the XY model in the 
rotated basis, we choose 

N N 

2 



Z i=l i>j 

In this case, the outlined approach gives for N > 5, f(N) ~ 0.25iV — 0.85. The simulation 
of the NN Ising model is achieved with 

m = ?E«f + E^°i. ( B14 ) 

Z i=l i>j 
N N-2 

A 2 = i E of - E J!,i + 2^t +2 , (bis) 

/ i=l i=l 

and, for > 5, f(N) ~ 0.015A^ - 0.035. 
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